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Abstract. Let _F be a real quadratic field with ring of integers and 
with class number 1. Let J" be a congruence subgroup of GL2(0). We 
describe a technique to compute the action of the Hecke operators on the 
cohomology H^{r\ C). For F real quadratic this cohomology group con- 
tains the cuspidal cohomology corresponding to cuspidal Hilbert modular 
forms of parallel weight 2. Hence this technique gives a way to compute 
the Hecke action on these Hilbert modular forms. 



1 Introduction 



1.1 



Let G be a reductive algebraic group defined over Q, and let F C G(Q) be an 
arithmetic subgroup. Let Y = r\X be the locally symmetric attached to G = 
G(M) and 7^, where X is the global symmetric space, and let be a local system 
on Y attached to a rational finite-dimensional complex representation of 7^. 
The cohomology H*{Y \Ai) plays an important role in number theory, through 
its connection with automorphic forms and (mostly conjectural) relationship to 
representations of the absolute Galois group Gal(Q/Q) (cf. [3,7,16,26]). This 
relationship is revealed in part through the action of the Hecke operators on the 
cohomology spaces. Hecke operators are endomorphisms induced from a family 
of correspondences associated to the pair (i^, G(Q)); the arithmetic nature of 
the cohomology is contained in the eigenvalues of these linear maps. 

For r C SL„ (Z) , modular symbols provide a concrete method to compute the 
Hecke eigenvalues in H'^{Y\ A^), where v = n{n+l) /2 — 1 is the top nonvanishing 
degree [8,23]. Using modular symbols many people have studied the arithmetic 
significance of this cohomology group, especially for n — 2 and 3 [3,6,7,12,25,26]; 
these are the only two values of n for which H'^{Y;AA) can contain cuspidal 
cohomology classes, in other words cohomology classes coming from cuspidal 
automorphic forms on GL(n). Another setting where automorphic cohomology 
has been profitably studied using modular symbols is that of C SL2(0), where 
is the ring of integers in a complex quadratic field [10, 11, 13, 22] . In this case Y 
is a three-dimensional hyperbolic orbifold; modular symbols allow investigation 
of H^{Y]M)^ which again contains cuspidal cohomology classes. 



1.2 

Now let F be a real quadratic field with ring of integers 0, and let G be the 
Q-group ReSi?/Q(GL2). Let F C G(Q) be a congruence subgroup. In this case we 



have X ~ x io x M, where Sj is the upper halfplane ( H2.2p . The locaUy symmetric 
space Y is topologically a circle bundle over a Hilbert modular surface, possibly 
with orbifold singularities if F has torsion. The cuspidal cohomology of Y is built 
from cuspidal Hilbert modular forms. Hence an algorithm to compute the Hecke 
eigenvalues on the cuspidal cohomology gives a topological technique to compute 
the Hecke eigenvalues of such forms. But in this case there is a big difference 
from the setting in ijl.ll the top degree cohomology occurs in degree v = A, 
but the cuspidal cohomology appears in degrees 2,30 Thus modular symbols 
cannot "see" the cuspidal Hilbert modular forms, and cannot directly be used 
to compute the Hecke eigenvalues. 

1.3 

In this article we discuss a technique, based on constructions in [17], that in prac- 
tice allows one to compute the Hecke action on the cohomology space H^{Y; C). 
Moreover it is easy to modify our technique to compute with other local sys- 
tems; all the geometric complexity occurs for trivial coefficients. Here we must 
stress the phrase in practice, since we cannot prove that our technique will actu- 
ally work. Nevertheless, the ideas in [17] have been successfully used in practice 
[4, 5] , and the modifications presented here have been extensively tested for 
F = Q(V2),Q(V3). 

The basic idea is the following. We first identify a finite topological model for 
H^{Y;C), the Voronoi reduced cocycles. This uses a generalization of Voronoi's 
reduction theory for positive definite quadratic forms [1,20], which constructs a 
/^-equivariant tessellation of X ( H2.2p . The Hecke operators do not act directly 
on this model, and to accommodate the Hecke translates of reduced cocycles we 
work with a larger model for the cohomology, the (infinite-dimensional) space 
Si{r) of 1-sharblies modulo F {^i^- The space Si{F) is part of a homological 
complex S^{F) with Hecke action that naturally computes the cohomology of F. 
Any Voronoi reduced cocycle in gives rise to a 1-sharbly cycle, which allows 
us to identify a finite dimensional subspace Sl'^'^{F) C Si{F). 

The main construction is then to take a general 1-sharbly cycle ^ and to 
modify it by subtracting an appropriate coboundary to obtain a homologous 
cycle 1^' that is closer to being Voronoi reduced (Q. By iterating this process, 
we eventually obtain a cycle that lies in our finite-dimensional subspace Sl^'^{F). 
Unfortunately, we are unable to prove that at each step the output cycle ^' is 
better than the input cycle ^, in other words that it is somehow "more reduced." 
However, in practice this always works. 

The passage from ^ to ^' is based on ideas found in [17], which describes 
an algorithm to compute the Hecke action on of congruence subgroups of 

^ The reader is probably more familiar with the case of G' — Res^/(jSL2. In this case 
the locally symmetric space is a Hilbert modular surface, and the cuspidal Hilbert 
modular forms contribute to H^. Our symmetric space is slightly larger since the 
real rank of G is larger than that of G'. However, regardless of whether one studies 
the Hilbert modular surface or our GL2 symmetric space, the cusp forms contribute 
to the cohomology in degree one below the top nonvanishing degree. 



SL4(Z). The common feature that this case has with that of subgroups of GL2(0) 
is that the cuspidal cohomology appears in the degree one less than the highest. 
This means that from our point of view the two cases are geometrically very 
similar. There are some complications, however, coming from the presence of 
non-torsion units in 0, complications leading to new phenomena requiring ideas 
not found in [17]. This is discussed in 21 We conclude the article by exhibiting 
the reduction of a 1-sharbly to a sum of Voronoi reduced 1-sharblies where the 
base field is Q(\/2) (©. 

We remark that there is another case sharing these same geometric features, 
namely that of subgroups of GL2{0k), where X is a complex quartic field. We 
are currently applying the algorithm in joint work with F. Hajir and D. Ramakr- 
ishnan for K = Q(C5) to compute the cohomology of congruence subgroups of 
GIj2{Ok) and to investigate the connections between automorphic cohomology 
and elliptic curves over K. Details of these cohomology computations, including 
some special features of the field K, will appear in [18]; the present paper focuses 
on the Hilbert modular case. 

Finally, we remark that there is a rather different method to compute the 
Hecke action on Hilbert modular forms using the Jacquet-Langlands correspon- 
dence. For details we refer to work of L. Dembele [14,15]. However, the Jacquet- 
Langlands technique works only with the complex cohomology of subgroups of 
GL2(0), whereas our method in principle allows one to compute with torsion 
classes in the cohomology. 

2 Background 
2.1 

Let F be a real quadratic field with class number 1. Let O C F denote the 
ring of integers. Let G be the Q-group ReSi?/Q(GL2) and let G — G(M) the 
corresponding group of real points. Let K d G he a, maximal compact subgroup, 
and let Ac be the identity component of the maximal Q-split torus in the center 
of G. Then the symmetric space associated to G is X = G/KAq- Let F C 
GL2(C) be a finite index subgroup. 

In §2.2l we present an explicit model of X in terms of positive-definite binary 
quadratic forms over F and construct a GL2(C')-equivariant tessellation of X 
following [1,20]. Section [2?3l recalls the sharbly complex [9,17,21]. 

2.2 Voronoi polyhedron 

Let ii, L2 be the two real embeddings of F into R. These maps give an isomor- 
phism F (g)Q M ~ , and more generally, an isomorphism 

G ^ GL2(M) X GL2(M). (1) 

When the meaning is clear from the context, we use ti, t2 to denote all such 
induced maps. In particular, ^ is the map 

9^M9),i2{9)). (2) 



Under this identification, Aq corresponds to {{rl,rl) \ r > 0}, where / is the 
2x2 identity matrix. 

Let C be the cone of real positive definite binary quadratic forms, viewed as 
a subset of V, the M- vector space of 2 x 2 real symmetric matrices. The usual 
action of GL2(M) on C is given by 

{g ■ (l)){v) = (l)Cgv), where g g GL2(R) and G C. (3) 

Equivalently, if is the symmetric matrix representing (jj, then g-cf) = gA^^g. In 
particular a coset gO{2) G GL2(M)/0(2) can be viewed as the positive definite 
quadratic form associated to the symmetric matrix g*^g. 

Let C = C X C. Then ^ and ^ define an action of G on C. Specifically, 
g ■ {4>i,4>2) = {ai,a2), where Ui is represented by Li{g)Acf,.ii{*g). Let (pQ denote 
the quadratic form represented by the identity matrix. Then the stabilizer in 
G of {4)o,(j)o) is a maximal compact subgroup K. The group Aq acts on C by 
positive real homotheties, and we have 

X = C/M>o = (C X C)/K>o -SjxS^xR, 

where is the upper halfplane. 

Let C denote the closure oiC inVxV. Each vector w e gives a rank 1 pos- 
itive semi-definite form w*w (here w is regarded as a column vector). Combined 
with Li and i2, we get a map L : — s- C given by 

Liv)^{L,iv)-XLM),'^2iv)-Xi2iv))). (4) 

Let R{v) be the ray R>o • L{v) C C. Note that 

L{cv) = iLi{cfLi{v),L2{cfL2{v)) 

so that if c e Q, then L{cv) G R{v), and in particular L{—v) — L{v). The set of 
rational boundary components Ci of C is the set of rays of the form R{v), v G F'^ 
[1]. These are the rays in C that correspond to the usual cusps of the Hilbert 
modular variety. 

Let A (ZV X V he the lattice 



A 



{(.1 (A),. 2(A)) 



A 



a c 
.cb. 



a, 6, cG O 



Then GL2(C') preserves A. 

Definition 1. The Voronoi polyhedron U is the closed convex hull in C of the 
points Ci n yl \ {0}. 

Since F has class number 1, one can show that any vertex of 7T has the form 
L{v) for V G 0^. We say that G 0^ is primitive if L{v) is a vertex of 77. Note 
that V is primitive only if L(v) is primitive in the usual sense as a lattice point 
in A. 



By construction GL2{0) acts on 77. By taking the cones on the faces of 77, 
one obtains a F -admissible decomposition of C for 7^ ~ Gh2{0) [1]. Essentially 
this means that the cones form a fan in C and that there are finitely many cones 
modulo the action of GL2{0). Since the action of GL2{0) commutes with the 
homotheties, this decomposition descends to a GL2(C')-equivariant tessellation 

oixE 

We call this decomposition the Voronoi decomposition. We call the cones 
defined by the faces of 77 Voronoi cones, and we refer to the cones corresponding 
to the facets of 77 as top cones. The sets tr H C, as cr ranges over all top cones, 
cover C. Given a point € C, there is a finite algorithm that computes which 
Voronoi cone contains [19]. 

For some explicit examples of the Voronoi decomposition over real quadratic 
fields, we refer to [24] (see also 3S])- 

2.3 The sharbly complex 

Let S'fc, fc > 0, be the 7^-module Ak/Ck, where Ak is the set of formal Z-linear 
sums of symbols [v] = [vi,--- ,i'fc+2], where each Vi is in F^, and Ck is the 
submodule generated by 

1- [va(l),--- ,'f<7(fc+2)] - sgn(cr)[ui, • • • ,Vk+2], 

2. [v,V2,-- - , Vk+2] - [w, W2, • • • Vk+2] if R{v) = 7?(w), and 

3. [v], if V is degenerate, i.e., if wi, • • • , are contained in a hyperplane. 

We define a boundary map d: Sk+i Sk by 

k+2 

d[vi,--- ,Vk+2] = ^(-l)'[^'l,- • • ,V^,■■■ ,Vk+2]- (5) 
i=l 

This makes S"* into a homological complex, called the sharbly complex [2]. 

The basis elements u = [vi,-- ■ , Ufe+2] are called k-sharblies. Notice that in 
our class number 1 setting, using the relations in Ck one can always find a 
representative for u with the Vi primitive. In particular, one can always arrange 
that each L{vi) is a vertex of 77. When such a representative is chosen, the Vi are 
unique up to multiplication by ±1. In this case the Vi — or by abuse of notation 
the L{vi) — are called the spanning vectors for u. 

Definition 2. A sharbly is Voronoi reduced if its spanning vectors are a subset 
of the vertices of a Voronoi cone. 

The geometric meaning of this notion is the following. Each sharbly u with 
spanning vectors Vi determines a closed cone cr(u) in C, by taking the cone 
generated by the points L{vi). Then u is reduced if and only if cr(u) is contained 
in some Voronoi cone. It is clear that there are finitely many Voronoi reduced 
sharblies modulo 7^. 

Using determinants, we can define a notion of size for 0-sharblies: 

^ If one applies this construction to F = Q, one obtains the Farey tessellation of 9), 
with tiles given by the SL2(Z)-orbit of the ideal geodesic triangle with vertices at 

0, 1, CX3. 



Definition 3. Given a 0-sharbly v, the size Size(v) of v is given by the absolute 
value of the norm determinant of the 2x2 matrix formed by spanning vectors 
for V. 

By construction Size takes values in Z>o. We remark that the size of a 0-sharbly 
V is related to whether or not v is Voronoi reduced, but that in general there 
exist Voronoi reduced 0-sharblies with size > 1 . 

The boundary map ^ commutes with the action of F, and we let S^{r) be 
the homological complex of coinvariants. Note that S^,{r) is infinitely generated 
as a Z_r-module. One can show 

Hkiis. ® c)(r)) ^ H^-'^ir; C) (6) 

(cf. [2]), with a similar result holding for cohomology with nontrivial coefficients. 
Moreover, there is a natural action of the Hecke operators on S^,{^) (cf. [17]). 
Thus to compute with iJ^(P;C), which will realize cuspidal Hilbert modular 
forms over F of weight (2,2), we work with 1-sharbly cycles. We note that 
the Voronoi reduced sharblies form a subcomplex of 5*(-r) finitely generated 
subcomplex that also computes the cohomology of -T as in © . This is our finite 
model for the cohomology of F. 

3 The reduction algorithm 

3.1 The strategy 

The general idea behind our algorithm is simple. To compute the action of a 
Hecke operator on the space of 1-sharbly cycles, it suffices to describe to an 
algorithm that writes a general 1-sharbly cycle as a sum of Voronoi reduced 1- 
sharblies. Now any basis 1-sharbly u contains three sub-O-sharblies (the edges of 
u), and the Voronoi reduced 1-sharblies tend to have edges of small size. Thus 
our first goal is to systematically replace all the 1-sharblies in a cycle with edges 
of large size with 1-sharblies having smaller size edges. This uses a variation 
of the classical modular symbol algorithm, although no continued fractions are 
involved. Eventually we produce a sum of 1-sharblies with all edges Voronoi 
reduced. However, having all three edges Voronoi reduced is (unfortunately) not 
a sufficient condition for a 1-sharbly to be Voronoi reduced]! Thus a different 
approach must be taken for such 1-sharblies to finally make the cycle Voronoi 
reduced. This is discussed further in Sj4l 

3.2 Lifts 

We begin by describing one technique to encode a 1-sharbly cycle using some 
mild extra data, namely that of a choice of lifts for its edges: 

^ This is quite different from what happens with classical modular symbols, and reflects 
the infinite units in 0. 



Definition 4 ([17]). ^2x2 matrix M with coefficients of F with columns Ai, 
A2 is said to be a lift of a 0-sharbly [u,v] if {R{Ai) , R{A2)} = {R{u),R{v)}. 

The idea behind the use of hfts is the foUowing. Suppose a linear combination 
of 1-sharblies ^ = ^a(u)u £ Si becomes a cycle in Si{r). Then its boundary 
must vanish modulo F. In the following algorithm, we attempt to pass from ^ 
to a "more reduced" sharbly ^' by modifying the edges of each u in the support 
of f . To guarantee that ^' is a cycle modulo F, we must make various choices 
in the course of the reduction F-equivariantly across the boundary of ^. This 
can be done by first choosing 2x2 integral matrices for each sub-O-sharbly of ^. 
We refer to [17] for more details and discussion. For the present exposition, we 
merely remark that we always view a 1-sharbly u = [vi,V2, V3] as a triangle with 
vertices labelled by the Vi and with a given (fixed) choice of lifts for each edge 
(Figured]). If two edges v, v' satisfy 7 • v = v', then we choose the corresponding 
lifts to satisfy 7M = M'. The point is that we can then work individually with 
1-sharblies enriched with lifts; we don't have to know explicitly the matrices in 
F that glue the 1-sharblies into a cycle modulo F. 

We emphasize that the lift matrices for any given 1-sharbly in the support 
of ^ are essentially forced on us by the requirement that f be a cycle modulo F. 
There is almost no flexibility in choosing them. Such matrices form an essential 
part of the input data for our algorithm. 




V2 Ml '"3. 



Fig. 1. A 1-sharbly with lifts 



3.3 Reducing points 

Definition 5. Let v be a 0-sharbly with spanning vectors {x,y}. Assume v is 
not Voronoi reduced. Then u G \ {0} is a reducing point for v if the following 
hold: 

L R{u)^R{x),R{y). 

2. L{u) is a vertex of the unique Voronoi cone a (not necessarily top-dimensional) 
containing the ray R{x + y). 

3. If X = ty for some t € F^ , then u is in the span of x. 

4- Of the vertices of a , the point u minimizes the sum of the sizes of the 0- 
sharblies [x^u\ and [u^y\. 



Given a non-Voronoi reduced 0-sharbly v — [x,y] and a reducing point u, we 
apply the relation 

[x,y] = [x,u] + [u,y] (7) 

in the hopes that the two new 0-sharblies created are closer to being Voronoi 
reduced. Note that choosing u uses the geometry of the Voronoi decomposition 
instead of (a variation of) the continued fraction algorithms of [8, 12,23]. Unfor- 
tunately we cannot guarantee that the new 0-sharblies on the right of ([7]) are 
better than v, but this is true in practice. 



3.4 P-invariance 

The reduction algorithm proceeds by picking reducing points for non- Voronoi 
reduced edges. We want to make sure that this is done -T-equivariantly; in other 
words that if two edges v, v' satisfy 7 • v = v', then if we choose u for v we want 
to make sure that we choose 7U for v'. 

We achieve this by making sure that the choice of reducing point for v only 
depends on the lift matrix M that labels v. The matrix is first put into normal 
form, which is a unique representative Mq of the coset GL2(0)\M. This is an 
analogue of Hermite normal form that incorporates the action of the units of 
O. There is a unique 0-sharbly associated to Mq; We choose a reducing point u 
for this 0-sharbly and translate it back to obtain a reducing point for v. Note 
that u need not be unique. However we can always make sure that the same 
u is chosen any time a given normal form Mq is encountered, for instance by 
choosing representatives of the Voronoi cones modulo GL2{0) and then fixing 
an ordering of their vertices. 

We now describe how Mq is constructed from M. Let fi^, be a fundamental 
domain for the action of {O^ ,■) on . For t € O, let i^+{t) be a fundamental 
domain for the action of {tO, +) on F. 

Definition 6. A nonzero matrix M e Mat2(i^) is in normal form if M has one 
of the following forms: 



1. 

2. 
3. 



ro b 

a b 

a b 
d 



, where b d f2^. 

, where a G /7» and b Cz F. 

, where a,d G and b G f2^(d). 



It is easy to check that the normal form for M is uniquely determined in the coset 

"a b' 



GLaCO) • M. To explicitly put M 



in normal form, the first step is to find 



. c d. 

7 G GL2(C') such that j-M is upper triangular. Such a 7 can be found after finite 
computations as follows. Let N : ^ R be defined by N{a) = | Normi?/Q(a)|. 
If 

< N{c) < N{a), 



then let a e O be an element of smallest distance from a/c. Let 





r 




1 -a 


7' = 


.10. 




.0 1 . 



Then 7' e GL2(0) and = 



with N{c') < N{c) and N{a') < N{a). 



Repeating this procedure will yield the desired result. 

After a reducing point is selected for v and the relation ([7]) is applied, we 
must choose lifts for the 0-sharblies on the right of ([7]). This we do as follows: 

Definition 7. Let [vi,V2] be a non-reduced 0-sharbly with lift matrix M and 
reducing point u. Then the inherited lift Mi for [vi, u] is the matrix obtained from 
M by keeping the column corresponding to Vi and replacing the other column by 
u. 



3.5 The algorithm 

Let T — [f 1, f2, W3] be a non-degenerate sharbly. Let Mi be the lifts of the edges 
of T as shown in Figure [TJ The method of subdividing the interior depends on 
the number of edges that are Voronoi reduced. After each subdivision, lift data 
is attached using inherited lifts for the exterior edges. The lift for each interior 
edge can be chosen arbitrarily as long as the same choice is made for the edge to 
which it is glued. We note that steps (I), (II), and (III.l) already appear in [17], 
but (ni.2) and (IV) are new subdivisions needed to deal with the complications 
of the units of O. 

(I) Three non-reduced edges. If none of the edges are Voronoi reduced, then 
we split each edge by choosing reducing points ui,U2, and M3. In addition, form 
three additional edges M2], ["2, M3], and [u3,mi]. We then replace T by the 
four 1-sharblies 

bl,1'2,W3] I > bl,'"3,W2] + [uz,V2,Ui] + [U2, "1,1^3] + [ui,1i2,U3]. 

(II) Two non-reduced edges. If only one edge is Voronoi reduced, then we 
split the other two edges by choosing reducing points ui and M3. We form two 
additional edges [^1,^3] and where £ is taken to be either or [w3,M3], 
whichever has smaller size. More precisely: 

1. If Size([i;i, ui]) < Size([w3, D3]), then we form two additional edges [1*1,^3] 
and [ui, ui], and replace T by the three 1-sharblies 

[vi,V2,Vz] I > [wi,M3,Ul] + [U'i,V2,Ui] + [wi , Ui , ^3] . 

2. Otherwise, we form two additional edges [ui, u^] and [^3, U3], and replace T 
by the three 1-sharblies 



(Ill) One non-reduced edge. If two edges are Voronoi reduced, then we split 
the other edge by choosing a reducing point ui. The next step depends on the 
configuration of {vi, W2, W3, ui}. 

1. If [v2, ui] or [wi, ^3] is not Voronoi reduced or V2 = tvi for some v G F, then 
we form one additional edge [vi , ui] and replace T by the two 1-sharblies 



2. Otherwise, a central point w is chosen. The central point w is chosen from 
the vertices of the top cone containing the barycenter of [vi, V2,V3, w] so that 
it maximizes the number of Voronoi reduced edges in the set 



We do not allow vi,V2 or V3 to be chosen as a central point. We form four 
additional edges w], [^2, w], [ui,w], and [^3, w] and replace T by the four 
1-sharblies 

[V1,V2,V3] I > [vi,V2,w] + [w,V2,Ui] + [w,Ui,V3] + [w,V3,Vi]. 

(IV) All edges Voronoi reduced. If all three edges are Voronoi reduced, but 
T is not Voronoi reduced, then a central point w is chosen. The central point w is 
chosen from the vertices of the top cone containing the barycenter of [vi , V2 , V3] 
so that it maximizes the sum + #P, where E is the set of Voronoi reduced 
edges in {[wi, w], [v2, w], [v3,w]} and P is the set of Voronoi reduced triangles in 
{[wi,W2,w], [v2,V3,w], [v3,vi,w]}. We do not allow vi,V2 or V3 to be chosen as 
a central point. We form three additional edges [wi,?/;], [v2,w], and [w3,w] and 
replace T by the three 1-sharblies 



4 Comments 

First, we emphasize that the reducing point u of Definition [S] works in practice 
to shrink the size of a 0-sharbly v, but we have no proof that it will do so. 
The difficulty is that Definition [5] chooses u using the geometry of the Voronoi 
polyhedron 77 and not the size of v directly. Moreover, our experience with 
examples shows that this use of the structure of 77 is essential to reduce the 
original 1-sharbly cycle (cf. §5.2|) . 

Next, as mentioned in §3.11 case (IV) is necessary: there are 1-sharblies T 
with all three edges Voronoi reduced, yet T is itself not Voronoi reduced. An 
example is given in the next section. The point is that in C the points L(v) and 
L{ev) are different if e is not a torsion unit, but after passing to the Hilbert 
modular surface L{v) and L(eu) define the same cusp. This means one can take 
a geodesic triangle A in the Hilbert modular surface with vertices at three cusps 




S = [v2,w], [v3,w], [ui,w]}. 



[V1,V2,V3] 



[vi,V2,w] + [w,V2,V3] + [w,V3,Vi]. 



that by any measure should be considered reduced, and can lift Zi to a 3-cone 
in the GL2-synimetric space that is far from being Voronoi reduced. 

Finally, the reduction algorithm can be viewed as a two stage process. When 
a 1-sharbly T has 2 or 3 non-reduced edges or 1 non-reduced edge and satisfies 
the criteria for case 1, then in some sense T is "far" from being Voronoi reduced. 
One tries to replace T by a sum of 1-sharblies that are more reduced in that 
the edges have smaller size. However, this process will not terminate in Voronoi 
reduced sharblies. In particular, if T is "close" to being Voronoi reduced, then 
one must use the geometry of the Voronoi cones more heavily. This is why we 
need the extra central point w in (III. 2) and (IV). 

For instance, suppose T = [wi , W2 , f 3] is a 1-sharbly with 1 non- reduced edge 
such that the criteria for (III. 2) are satisfied when the reducing point is chosen. 
One can view choosing the central point and doing the additional split as first 
moving the bad edge to the interior of the triangle, where the splitting no longer 
needs to be r'-invariant. The additional freedom allows one to make a better 
choice. Indeed, without the central point chosen wisely, this does lead to some 
problems. In particular, there are examples where [ui, Ui] is not Voronoi reduced, 
and the choice of the reducing point for this edge is V2, leading to a repeating 
behavior. Thus the distinction had to be made. 

5 The case F = Q(a/2) 
5.1 

Let F = Q(-\/2) and let e = 1 + v^, a fundamental unit of norm —1. Computa- 
tions of H. Ong [24, Theorem 4.1.1] with positive definite binary quadratic forms 
over F allow us to describe the Voronoi polyhedron 77 and thus the Voronoi de- 
composition of C: 

Proposition 1 ([24, Theorem 4.1.1]). Modulo the action o/GL2(0), there 
are two inequivalent top Voronoi cones. The corresponding facets of 11 have 6 
and 12 vertices, respectively. 

We fix once and for all representative 6-dimensional cones ^0 and Ai. To 
describe these cones, we give sets of points S* C such that the points {L{v) \ 
V & S} are the vertices of the corresponding face of 71. Let ei, 62 be the canonical 
basis of O"^. Then we can take ^0 to correspond to the 6 points 

ei, 62, ei - 62, eei, £62, £(61 - 62), 

and Aq to correspond to the 12 points 

61, 62, £61, £62, 61 — 62, 61 -I- £62, 62 -I- £61, £(61 + 62), a, (3, ea,eP, 

where a = ei — \/2e2, /? = 62 — V^ei. Since Ai is not a simplicial cone, there 
exist basis sharblies that are Voronoi reduced but do not correspond to Voronoi 
cones. 



Now we consider cones of lower dimension. Modulo GL2{0), every 2-dimensional 
Voronoi cone either lies in C \ C or is equivalent to the cone corresponding to 
{61,62}. The GL2(0)-orbits of 3-dimensional Voronoi cones are represented by 
{61, 62} U J7, where U ranges over 

{ei - 62}, {eei}, {e(ei - 62)}, {ei - V2e2, 62 - V2ei}, {ei + £62}. 

Note that all but one of the 3-cones are simplicial. 

5.2 

Now we consider reducing a 1-sharbly T. Let us represent T by a 2 x 3 matrix 
whose columns are the spanning vectors of T. We take T to be 

r\/2 + 34V2 + 4 3\/2-4 " 
[ V2 5V2-I -3V2-5J' 

and we choose arbitrary initial lifts for the edges of T. This data is typical of 
what one encounters when trying to reduce a 1-sharbly cycle modulo F. 

The input 1-sharbly T has 3 non-reduced edges with edge sizes given by the 
vector [5299, 529, 199]. The first pass of the algorithm follows (I) and splits all 3 
edges, replacing T by the sum Si + S2 + S3 + S4, where 

4V24-4 -V2-1' 

5V2-1 -V2-1 -V2 y 

1-V2-1' 
-V2-IO -V2 ■ 

We compute that Size(5i) = [2,2,8], Size(52) = [1,1,16], Size(5'3) = [1,2,7], 
and Size(S'4) = [2, 1, 1]. Notice that the algorithm replaces T by a sum of shar- 
blies with edges of significantly smaller size. This kind of performance is typical, 
and looks similar to the performance of the usual continued fraction algorithm 
over Z. Note also that S4, which is the 1-sharbly spanned by the three reduc- 
ing points of the edges T, also has edges of very small size. This reflects our 
use of Definition [5] to choose the reducing points; choosing them without using 
the geometry of U often leads to bad performance in the construction of this 
1-sharbly. 

Now 5*4 has 3 Voronoi reduced edges, but is itself not Voronoi reduced. The 
algorithm follows (IV), replaces S4 by Ri + R2 + R3, and now each Ri is Voronoi 
reduced. 

The remaining 1-sharblies Si, S2, and 53 have only 1 non-reduced edge. 
They are almost reduced in the sense that they satisfy the criteria for (III. 2). 
The algorithm replaces by Oi-f O2 + O3-I-O4, where Oi and O2 are degenerate 
and O3 and O4 are Voronoi reduced. The 1-sharbly S2 is replaced by a Pi -I- P2 -|- 
P3 + P4, and each Pi is Voronoi reduced. S3 is replaced hy Qi + Q2 + Q3 + Q4, 



Si 



V2 + 3 
V2 



V2- 1 1 
-^/2 



^2 



S3 = 



3\/2-4 1 
-3\/2 - 5 -V2 - 1 



^4 = 



where Qi and Q2 are degenerate, Q3 is Voronol reduced, and Q4 is not Voronoi 
reduced. This 1-sharbly is given by 

^ _ r-\/2 + l 3^/2-4 " 
~ [ 2^/2 + 3 -V2 - 1 -3V2 - 5 

and has 3 Voronoi reduced edges. Once again the algorithm is in case (IV), and 

replaces Q4 by a sum A^^i + + ^3 of Voronoi reduced sharblics. 

To summarize, the final output of the reduction algorithm applied to T is a 
sum 

Ni+ N2 + N3 + O3 + O4 + Pi+ P2 + P3 + P4 + Q3 + Ri+ R2 + R3, where 



^ I + 1 

^ I 2V2 + 3 -V2-I -2V2- 



N2 



N3 



O4 



Pi 



P2 



Pa = 



3V2-4 -V2 + 1 

-3\/2 - 5 2^/2 + 3 -2\/2 - 3 

-\/2- 1 1 + 5 

-1 V2 

2V2 + 3 1 

\/2 + 2 -2\/2 + 2 -\/2- 1 

2\/2 + 3 -\^- 1 4\/2 + 4 
\/2 + 2 -1/2 5^2-1 

\/2 + l 



3\/2-4 
V2 - 1 -3V2 - 5 -2V2 - 3 

-72- 1 -v^- 1 1 



O3 



-1 



-\/2 

2V2 + 34V2 + 4 1 
V2 + 2 5V2- 1 -2V2 + 2 

2V2 + 3 -V2-1 

\/2 + 2 -V2-I -\/2 

-\/2 + ll 
2a/2 + 3 -V2-I 



Q 



Ri = 



R2 = 



R3 




2^/2 + 3 -V2 - 1 -2^/2 - 3 

3V2-4 

\/2 - 1 -3\/2 - 5 -2V2 - 3 

3^2-4 -a/2 + 1 
-3^2-5 2a/2 + 3 -2\/2-3 



and 



and each of the above is Voronoi reduced. Some of these 1-sharblics correspond 
to Voronoi cones and some don't. In particular, one can check that the spanning 
vectors for P3, P4, Ri, and A^i do form Voronoi cones, and all others don't. 
However, the spanning vectors of O3 and O4 almost do, in the sense that they 
are subsets of 3-dimensional Voronoi cones with four vertices. 



References 



1. A. Ash, Deformation retracts with lowest possible dimension of arithmetic quotients 
of self-adjoint homogeneous cones, Math. Ann. 225 (1977), 69-76. 



2. , Unstable cohomology ofSL{n,0), J. Algebra 167 (1994), no. 2, 330-342. 

3. A. Ash, D. Grayson, and P. Green, Computations of cuspidal cohomology of con- 
gruence subgroups of SL-i{7i), J. Number Theory 19 (1984), 412-436. 

4. A. Ash, P. E. Gunnells, and M. McGonnell, Cohomology of congruence subgroups 
of SL4{Z) 11, submitted to J. Number Theory. 

5. , Cohomology of congruence subgroups of SL4{Z), J. Number Theory 94 

(2002), 181-212. 

6. A. Ash and M. McConnell, Experimental indications of three-dimensional Galois 
representations from the cohomology of SL{3,Z), Experiment. Math. 1 (1992), 
no. 3, 209-223. ^ 

7. A. Ash, R. Pinch, and R. Taylor, An Aa extension of Q attached to a non-selfdual 
automorphic form on GL{3), Math. Ann. 291 (1991), 753-766. 

8. A. Ash and L. Rudolph, The modular symbol and continued fractions in higher 
dimensions, Invent. Math. 55 (1979), 241-250. 

9. A. Ash, Unstable cohomology ofSL{n,0), J. Algebra 167 (1994), no. 2, 330-342. 

10. J. Bygott, Modular forms and modular symbols over imaginary quadratic fields, 
Ph.D. thesis, Exeter, 1999. 

11. J. E. Cremona, Hyperbolic tessellations, modular symbols, and elliptic curves over 
complex quadratic fields, Compositio Math. 51 (1984), no. 3, 275 324. 

12. , Algorithms for modular elliptic curves, second cd., Cambridge University 

Press, Cambridge, 1997. 

13. J. E. Cremona and E. Whitley, Periods of cusp forms and elliptic curves over 
imaginary quadratic fields. Math. Comp. 62 (1994), no. 205, 407-429. 

14. L. Dembelc, Explicit computations of Hilbert modular forms on Q(-s/5), Experi- 
ment. Math. 14 (2005), no. 4, 457-466. 

15. , Quaternionic Manin symbols, Brandt matrices, and Hilbert modular forms. 

Math. Comp. 76 (2007), no. 258, 1039-1057 (electronic). 

16. J. Pranke, Harmonic analysis in weighted L2-spaces, Ann. Sci. Ecole Norm. Sup. 
(4) 31 (1998), no. 2, 181-279. 

17. P. E. Gunnells, Computing Hecke eigenvalues below the cohomological dimension. 
Experiment. Math. 9 (2000), no. 3, 351 367. 

18. P. E. Gunnells and D. Yasaki, Computing Hecke operators on modular forms over 
real quadratic and complex quartic fields, in preparation. 

19. P. E. Gunnells, Modular symbols for Q-rank one groups and Voronoi reduction, J. 
Number Theory 75 (1999), no. 2, 198-219. 

20. M. Koccher, Beitrdge zu einer Reduktionstheorie in Positivitatsbereichen. I, Math. 
Ann. 141 (1960), 384-432. 

21. R. Lee and R. H. Szczarba, On the homology and cohomology of congruence sub- 
groups, Invent. Math. 33 (1976), no. 1, 15-53. 

22. M. Lingham, Modular forms and elliptic curves over imaginary quadratic fields, 
Ph.D. thesis, Nottingham, 2005. 

23. Y.-I. Manin, Parabolic points and zeta-functions of modular curves. Math. USSR 
Izvestija 6 (1972), no. 1, 19-63. 

24. H. E. Ong, Perfect quadratic forms over real- quadratic number fields, Geom. Ded- 
icata 20 (1986), no. 1, 51-77. 

25. W. Stein, Modular forms, a computational approach. Graduate Studies in Math- 
ematics, vol. 79, American Mathematical Society, Providence, RI, 2007, With an 
appendix by P. E. Gunnells. 

26. B. van Geemen, W. van der Kallen, J. Top, and A. Verberkmoes, Hecke eigen- 
forms in the cohomology of congruence subgroups of SL(3,Z), Experiment. Math. 
6 (1997), no. 2, 163-174. 



